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Abstract 

We consider the dynamics of a field coupled to a harmonic crystal with n compo- 
nents in dimension d, d, n > 1 . The crystal and the dynamics are translation-invariant 
with respect to the subgroup 7L d of IR d . The initial data is a random function with 
a finite mean density of energy which also satisfies a Rosenblatt- or Ibragimov-Linnik- 
type mixing condition. Moreover, initial correlation functions are translation-invariant 
with respect to the discrete subgroup 7L d . We study the distribution fx t of the solution 
at time t G M . The main result is the convergence of fit to a Gaussian measure as 
t — > 00 , where /Uqo is translation-invariant with respect to the subgroup 2 d . 
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1 Introduction 



The paper concerns problems of long-time convergence to an equilibrium distribution in 
a coupled system which is similar to the Born-Oppenheimer model of a solid state. In 
[5, 6, 7, 10] we have started the convergence analysis for partial differential equations of 
hyperbolic type in IR d . In [8, 9] we have extended the results to harmonic crystals. 

Here we treat a harmonic crystal coupled to a scalar Klein-Gordon field. In this case, 
the corresponding problem in the unit cell is an infinite-dimensional Schrodinger operator, 
whereas in [8, 9] (and in [5, 6, 7, 10]) it was a finite-dimensional matrix. This situation 
usually arises in the solid-state problems similar to that for the Schrodinger equation with 
space-periodic potential [18]. The main novelty in our methods consists in that they yield 
exact estimates of trace norms for the problem in the unit cell. 

We assume that an initial state Y of the coupled system is a random element of a 
Hilbert phase space £ , see Definition 2.4. The distribution of Yq is a probability measure 
fiQ of mean zero satisfying conditions S1-S3. In particular, the measure fiQ is invariant 
with respect to translations by vectors of Z d . For a given t e IR , we denote by [i t the 
probability measure defining the distribution of the solution Y(t) to the dynamical equations 
with random initial state Y . We study the asymptotics of ji t as t — > ±oo . 

Our main result gives the (weak) convergence of the measures ji t to a limit measure ji^ , 



The measure //oo is Gaussian and translation-invariant with respect to the group Z . We 
give the explicit formulas for the covariance of the measure //oo . The dynamical group is 
ergodic and mixing with respect to the limit measure //qo . Similar results hold as t — > — oo 
because the dynamics is time-reversible. 

Similar results have been established in [1, 24] for one-dimensional chains of harmonic 
oscillators (with d — 1) and in [11, 13, 16, 21] for one-dimensional chains of anharmonic 
oscillators coupled to heat baths. For d -dimensional harmonic crystals, with d > 1 , the 
convergence (1.1) was proved in [8, 9, 17]. The mixing condition was first introduced by 
R. Dobrushin and Yu. Suhov for an ideal gas in [3]. The condition can replace the (quasi-) 
ergodic hypothesis when proving the convergence to the equilibrium distribution, and this 
plays a crucial role in our approach. Developing a Bernstein-type approach, we have proved 
the convergence for the wave and Klein-Gordon equations and for harmonic crystals with 
translation-invariant initial measures, [5, 6, 8]. In [7, 9, 10] we have extended the results 
to two-temperature initial measures. The present paper extends our previous results to the 
scalar Klein-Gordon field coupled to the nearest neighbor crystal. 

Let us outline our main result and the strategy of the proof. (For the formal definitions 
and statements, see Section 2.) Consider the Hamilton system with the following Hamilto- 
nian functional: 



t 



OO. 



(1.1) 




+ III {Y.Hk + e 3 )-u(k)f + \v(k)f + vZ\u(k)\ 2 ) 



+ 



kez d 




(1.2) 
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involving a real scalar field ip(x) and its momentum n(x) , x G IR , coupled to a "simple 
lattice" described by the deviations u(k) G IR n of the "atoms" and their velocities v(k) G 
IR™, k G 7Z. d . The symbol R{x) stands for a IR n -valued function and G Z d for the 
vector with the coordinates e*- := 5*- . Taking the variational derivatives of if(^, w, 7r, v) , we 
formally obtain the following system for x G M d and fc G 7L d : 





5H 

5n 


= n(x,t), 




it(k,t) = 


dH 


= v(k,t), 




dv 




•k(x,t) = 


SH 

!hp 


= (A-m 2 )i/j(x,t) 


- E u{k',t)-R(x-k'), 
k'ez d 


v(h,t) = 


dH 

du 


= (A L -v*)u(k,t) 


-J R(x' - k)ip(x',t) dx'. 



(1.3) 



Here mo, > , and denotes for the discrete Laplace operator on the lattice 2 d , 

A L u{k) := ]T (u(k + e)-u(k)). 

e,\e\=l 

Note that for n = d and R{x) = — Vp(x) , the interaction term in the Hamiltonian is the 
linearized Pauli-Fierz approximation of the translation-invariant coupling 

J p(x - k - u(h))ip(x) dx. (1.4) 

k 

A similar model was analyzed by Born and Oppenheimer [2] as a model of a solid state (cou- 
pled Maxwell-Schrodinger equations for electrons in the harmonic crystal; see, e.g., [18]). 
The traditional analysis of the coupled field-crystal system (1.3) is based on an iterative per- 
turbation procedure using the adiabatic approximation. Namely, in the zero approximation, 
the crystal and the (electron) field are discoupled. In the first one, the electron field defines 
a slow displacement of nuclei. The displacements give the corresponding contribution to 
the field via the static Coulombic potentials, which means a non-relativistic approximation, 
etc. The iterations converge if the motion of the nuclei is sufficiently slow, i.e., the nuclei are 
rather heavy as compared with the electrons. A similar procedure applies to the corresponding 
stationary problem of finding the dispersion relations. 

Our analysis of the dispersion relations is a bit different and holds for small displacements. 
Namely, we linearize the translation-invariant coupling (1.4) at the zero displacements of 
the nuclei and obtain the equations (1.3) corresponding to the Pauli-Fierz approximation. 
On the other hand, we analyze the dispersion relations of the linearized equations without 
any adiabatic or non-relativistic approximation. We give an exact nonperturbative spectral 
analysis of the coupled system (1.3). 

We study the Cauchy problem for the system (1.3) with the initial data 

ip(x,0) =VoOe), 7r(x,0) = 7r (x), x G IR d , , s 

u(k,0) =uo(k), v(k,0) =v (k), ke7L d . 1 } 
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Let us write 

ip° := if), ip l := it, u° := u, u 1 := v, 

vm ._ (Y o m Y i (f ss Y°(t) := ^°(x,t),u°(k,t)) := (ip(x, t), u(k, t)), 

i [i). v w,i Y l {t):={^ l {x,t),u 1 {k,t)):={7i{x,t),v{k,t)). 1 J 

In other words, Y(-, t) are functions defined on the disjoint union IP := IR d U Z d , 

In this case, the system (1.3), (1.5) becomes a dynamical problem of the form 

Y(t) = AY(t), t e IR; Y(0) = Y . (1.7) 
Here Y = (if) ,u ,7i ,v ) and 

where Su(x) = J2kez d R( x ~ k)u(k) , S*ip(k) = / R(x — k)ip(x) dx , and 

(^,^) L2(Rd) = (5>, M ) p(zd)]n , ^eL 2 (IR d ), M G[/ 2 (Z d )]". 
We assume that the initial datum Y is a random function, and the initial correlation matrix 

Qo(p,p') := e(y ( P )®Y ( P ')) : p,p' e P, 

is translation invariant with respect to translations by Z d , i.e., 

Q ( p + k,p' + k) = Q (p,p'), p,p'eP, (1.9) 

for any A; G Z d . We also assume that the initial mean energy densities are uniformly 
bounded, 

e F (x) := E(\Vip (x)\ 2 + \M%)\ 2 + Ko(^)| 2 ) < e F < oo, a. a. a; G IR ', (1.10) 
e L := E(\u (k)\ 2 + \v (k)\ 2 )<oo, ke7L d . (1.11) 

Finally, we assume that the measure /io satisfies a mixing condition of a Rosenblatt- or 
Ibragimov-Linnik type, which means that 

Y (p) and Y (p') are asymptotically independent as \p — p'| — > oo. (1-12) 

Our main result gives the (weak) convergence (1.1) of \i t to a limit measure fi^ , which is 
a stationary Gaussian probability measure. 

Let us comment on the methods of the proof. The key role in our proof is played by 
the standard reduction of system (1.7) to the Bloch problem on the torus. Namely, we split 
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x G IR d in the form x = k + y , k G Z d , y G i^i := [0, l] d , and apply the Fourier transform 
Fk-^e to the solution Y{k, t) :— \^{k + y, t), u(k, t), n(k + y, t), v(k, t) ) , 



Y(9, t) := F k ^ e Y(k, f) = £ e ife V (fc, t) = (^(0, </, i), u(9, t), n(9, y, t),v(9, t)), 9e IR d , 

kez d 

which is a version of the Bloch-Floquet transform. The functions ip , ft are periodic with 
respect to 9 and quasi-periodic with respect to y , i.e., 

#(e,y + m,t) = e- im9 ilj(6,y,t), ft(9, y + m, t) = e~ im6 ^{9, y, t), m G H d . 

Further, introduce the Zak transform of Y(-, t) (which is also known as Lifshitz-Gelfand-Zak 
transform) (cf [18, p. 5]) as 

ZY(-,t) = Y n (9,t) := (MO,y,t),u(e,t),7t n (e,y,t),v(0,t)), (1.13) 

where ipu(9,y,t) := e iyd ip(9,y,t) and ftu(9,y,t) := e ty6 ft(9,y,t) are periodic functions with 
respect to y (and quasi-periodic with respect to 9). Denote by Tf := IR d /Z d the real 
d -torus and write 1Z := Tf U {0} . Set 

Y„(i,r,t) = Yn(i,t) := { <^ V 'm rtf B ' ' )) ' 'T ^ 

1 (u(0,t),u(0,t)), r = 0. 

Problem (1.7) is now equivalent to the problem on the unit torus y G Tf with the parameter 
9eK d = [0, 27r] d , 



y n (0,f) =.A(0)Yh(0,t), t G H 

Vh(0, 0)=F nW 



9eK d . (1.14) 



Here 

= ( _ J ) , (1-15) 
and 7Y(0) := ZHZ^ 1 is the "Schrodinger operator" on the torus Tf , 

where 

ul{9) :=2(l-cos6* 1 ) + ... + 2(l-cos6) d ) + z/ 2 , (1.17) 
(S(6)u(-j)(e,y):=R u (e,y)-u(e), (S*(9)M9, -))(9) := [ ^(-6^)^(6, y) dy, (1.18) 

<tfn(0, •), (5(fl)fi)(fl, •) W*) = (S* Wn)(0) • «(0), tfn(0, •) G tf 1 ^), 5(0) G (D n . 
Then, formally, 

Yn(9,t) = e A ^% n (9), 9 G K d . (1.19) 

To justify the definition of the exponential, we note that 7i.(0) is a self-adjoint operator with 
a discrete spectrum. Indeed, if R — , then this follows from elliptic theory, and, if R ^ , 
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then the operators S(9) and S*(9) are finite-dimensional for a fixed 9. We assume that 
Ti(9) > (condition R2) which corresponds to the hyperbolicity of problem (1.3). 

Note that in [8, 9], we considered the harmonic crystal without any field. In this case, 
the operator A(9) is a finite-dimensional matrix. 

Let us prove the convergence (1.1) by using the strategy of [5]-[10] in the following three 
steps. 

I. The family of measures /i t , t > , is weakly compact in an appropriate Frechet space. 

II. The correlation functions converge to a limit, 

Qt(p,p') ■= J (Y(p) Y(p')) fr(dY) - Qocfep'), t - oo, p,p'e IP. (1.20) 

III. The characteristic functionals converge to a Gaussian functional, 

tk{Z):= Je i ^ fJh (dY)^exp{-^Q 00 {Z,Z)}, t- oo, (1.21) 

where Z is an arbitrary element of the dual space and is a quadratic form. 

Property I follows from the Prokhorov Theorem. First, let us prove the uniform bound 
(2.15) for the mean local energy in fi t ■ To this end, we shall show that the operator 
(p, l q t j (9)n j ^j . _ Q i is of trace class, where q t j (9) represents the covariance of the measure /j, t 

in the Zak transform (see (3.7)) and Q = Q(9) := \fTt{&) . Moreover, we derive the uniform 
bound 

sup sup trfft^fl)^') < oo. (1.22) 

t>o ee[o,2n] d 

This implies the compactness of fit by the Prokhorov theorem (when applying Sobolev's 
embedding theorem as in [5]). 

To derive property II, we study oscillatory integrals in the Zak transform by developing 
our cutting strategy intoduced in [8]. Namely, we rewrite (1.20) in the form 

Q t (Z,Z)^Q 00 (Z,Z), t^oo, (1.23) 

where Qt(Z, Z) stands for the correlation quadratic form for the measure /i t . Further, we 
prove formula (1.23) for Z G V° as follows: by the definition of T>° , the Zak transform 
Zu{9) vanishes in a neighborhood of a "critical set" C C K d . In particular, the set C 
includes all points 9 G K d with a degenerate Hessian of uji(9) and the points for which the 
function uoi(9) is non-smooth. One can cut off the critical set C by the following two crucial 
observations: (i) mes C = and (ii) the initial correlation quadratic form is continuous in L 2 
due to the mixing condition. The continuity follows from the spatial decay of the correlation 
functions in accordance with the well-known Shur lemma. 

Similarly, we first prove property III for Z <E V° and then extend it to all Z G T> . For 
Z G T>° , we use a version of the S.N. Bernstein "room-corridor" technique (cf. [5, 8]). This 
leads to a representation of the solution as the sum of weakly dependent random variables. 
Then (1.21) follows from the Central Limit Theorem under a Lindeberg-type condition. 

Let us comment on the two main technical novelties of our paper. The first of them is the 
bound (1.22), which then ensures compactness. We derive formula (1.22) in Section 4 directly 
from our assumption concerning the finiteness of the mean energy density (1.10), (1.11). 
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The derivation uses the technique of trace class operators [23], which enables us to avoid 
additional continuity conditions for higher-order derivatives of the correlation functions. An 
essential ingredient of the proof is the "unitary trick" (4.6), which is a natural consequence 
of the Hamiltonian structure of system (1.3). The second main novelty is the bound (2.2) for 
the dynamics in weighted norms. In the Zak transform, the weighted norms become Sobolev 
norms with negative index. We derive (2.2) in Appendix A, by using duality arguments, 
from the corresponding bounds for the derivatives of the exponential (1.19). The bounds for 
the derivatives follow by differentiating the dynamical equations. 

Let us comment on our conditions El and E2. The conditions are natural generalizations 
of similar conditions in [8, 9]. Condition El enables us to apply the stationary phase method 
to the oscillatory integral representation for the covariance. It provides that the stationary 
points of the phase functions are non-degenerate. 

The paper is organized as follows. In Section 2 we formally state our main result. The 
compactness (Property I) is established in Section 4, the convergence (1.20) in Section 6, and 
the convergence (1.21) in Section 7. In Section 8, mixing properties for the limit measures 
are proved. Appendix A concerns the dynamics in the Fourier transform, in Appendix B we 
analyze the crossing points of the dispersion relations, and in Appendix C we discuss the 
covariance in the spectral representation. 

2 Main results 
2.1 Notation 

We assume that the initial data Y are given by an element of the real phase space 8 defined 
below. 

Definition 2.1 Let H s ' a = H s ' a (M d ) , s G M , a G M , be the Hilbert space of distributions 
ip G S"(IR d ) with finite norm 



IMka = ll(^) a A>|| L2(]Rd) < 00. 

For i) G D = C °°(IR d ) ; write Fip(£) = f e* x tj)(x)dx. Let A s iP := i^((£)W)) and 



(x) := J\x\ 2 + 1 , where i[) := Ftp stands for the Fourier transform of a tempered distribution 



Remark 2.2 For s = 0, 1, 2, ... , the space H s ' a (JR d ) is the Hilbert space of real-valued 
functions ip(x) with finite norm 



which is equivalent to H^llsa ■ 

Definition 2.3 Let L a , a G JR , be the Hilbert space of vector functions u(k) G IR™ , 





k G Z d , with finite norm 




2 _ 



(l + |£f)X£;)| 2 < oo. 



a 



kez d 
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Definition 2.4 Let £ s > a := H 1+s ' a (JR d )®L a ®H s ' a (]R d )®L a be the Hilbert space of vectors 
Y = {jp, u, n, v) with finite norm 

lylla^U\\l + s,a+\Hl+\\A\la+\Hl- 

Choose some a , a < —d/2 . Assume that Y G £ := £ 0,a . 

Using the standard technique of pseudo-differential operators and Sobolev's Theorem 
(see, e.g., [14]), one can prove that £°' a = £ C £ s,lS for every s < and (5 < a , and the 
embedding is compact. 

Definition 2.5 TTie phase space of problem (1.7) is £ := £°' a , a < —d/2 . 

Introduce the space H[ := H S (T?) ©<D n , s G IR, where H s (Tf) stands for the Sobolev 
space. 

We assume that the following conditions hold for the real- valued coupling vector function 
R(x): 

Rl. ReC°°(JR d ) and \R(x)\ < R exp(— e\x\) with some e > and some R < oo . 

R2. The operator 7Y(0) is positive definite for 9 G if d = [0,27r] d . This is equivalent to the 

uniform bound 

(X°,7i:(#)X ) > k 2 \\XXh\ for ^° e itf, 9 G (2.1) 
where k > is a constant and (•, •) stands for the inner product in H® (see (3.13)). 

Remark 2.6 i) Condition R2 ensures that the operator iA(9) is self-adjoint with respect 
to the energy inner product. This corresponds to the hyperbolicity of problem (1.3). 

ii) Condition R2 holds, in particular, if the following condition R2' holds (see Remark 
9.3): 

R2'. / I E R(k + y)\ 2 dy<u 2 ml/2. 
J[0 ' 1] fcez d 

ra) Condition R2' holds for functions R satisfying condition Rl wrat/j <C 1 . 

Proposition 2.7 Lei conditions Rl and R2 /ioW. T/ien (i) for any Y G £ , there exists a 

unique solution Y(t) G C(JR,£) to the Cauchy problem (1.7). 

(ii) The operator W(t) : Y \— > F(t) zs continuous in £ for any t G IR, 

sup |||W(t)yb|||o,« < C(T)|||Yb||| 0> a (2.2) 

|t|<T 

i/ ck is even and a < — 2 . 

Proof. fi) Loca/ existence. Introduce the matrices 

■^^(-Wo o)' Wo= (~ A ()" m ° -A L ° +I / 2 )' W) 

Then problem (1.3) can be rewritten as the Duhamel integral 

t 

Y{t) = e A "% + J e Ao(t - s) BY(s)ds, 
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where 



BY = (o, 0, - ]T u(k')R(x - k'), - [ R(x' - k)iP(x') dx') , Y = (ip, u, vr, v). 



k'ez d 



Condition Rl implies that |||Sy(s)||| , Q < C|||r(s)||| , Q . Further, for < s < t < T , we 
obtain 

le^^BY(s)h^<C(T)ms)l^ 

(see, e.g., [8]). Hence, 

max \\\Y(t)\l , a < C(T)|||y ||o,« + TC(T) max ||y( a )||o,« < (T + l)C(T) max \\Y{s)\\,, a . 

We choose a T > so that (T + l)C(T) < 1. Then the contraction mapping principle 
implies the existence of a unique solution Y{t) e C([0,T];5) . T/ie global existence follows 
from the bound (2.2). 

(ii) The bounds (2.2) are proved in Corollary 9.6. ■ 



Conditions Rl and R2 imply that, for a fixed 9 G K d , the operator H(9) is positive 
definite and self-adjoint in H® and its spectrum is discrete. Introduce the Hermitian positive- 
definite operator 

£1(6) := yfR{&) > 0. 

Denote by u>i(9) > and Fi(9,-), I = 1,2,..., the eigenvalues ("Bloch bands") and the 
orthonormal eigenvectors ("Bloch functions") of the operator Q(9) in H ° , respectively. 
Note that F[(9, •) e iff := C°°(Tf) © C n , because these are eigenfunctions of the elliptic 
operator 7Y(#) . 

As is well known, the functions uoi(-) and Fi{-,r) are real-analytic outside the set of the 
"crossing" points , where u)i{9*) = ui'(6#) for some I ^ V . However, the functions are not 
smooth at the crossing points in general if wi(6) ^ wi>(9) . Therefore, we need the following 
lemma, which is proved in Appendix B. 

Lemma 2.8 (cf. [26]) There exists a closed subset C* C K d such that (i) the Lebesgue 
measure of C* vanishes, 

mesC* = 0. (2.4) 

(ii) For every point O G K d \ and N G IN , there exists a neighborhood 0(0) C K d \ 
such that each of the functions uii{9) and Fi{9,-) , I — 1, . . . , N , can be chosen to be real- 
analytic on 0(Q) . 

(Hi) The eigenvalues u>i(9) have constant multiplicity in 0(Q) , i.e., one can enumerate 
them in such a way that 

u 1 (e) = ..=u ri (6), iv ri+1 (9) = ...=u; r2 (9),..., (2.5) 
u r M^rM if a^v, r a ,r v >\, (2.6) 

for any 9 G (9(6) . 
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Corollary 2.9 The spectral decomposition holds, 

+00 

= 5>,(0)W), ^o(e), (2.7) 
1=1 

where Pi (9) are the orthogonal projectors in onto the linear span of Fi{9, ■) , and Pi{9) 
and u>i(9) depend on 9 G 0(Q) analytically. 

Assume that system (1.7) satisfies the following conditions El and E2. For every G 
K d \ C* : 

El £>,(0) ^ 0, I = 1,2,..., where D,(0) := det(|^)*. =i , 9 G (9(6) , and 0(6) is 
defined in Lemma 2.8. 
Write 

C,:= (J {0gO(9): A(0) = O}, Z = 1,2,.... 
The following lemma is also proved in Appendix B. 

Lemma 2.10 Let conditions Rl and R2 hold. Then mes Ci = , 1 — 1,2, 

E2 For each I ^ I 1 , the identity u>i(9) — u>i>(9) = const_ , 9 G 0{Q) , cannot hold for any 
constant const _ 7^ , and the identity uji{9) + uji>{9) = const + cannot hold for any constant 
const + 7^ . 

Condition E2 could be considerably weakened (cf. [8, Remark 2.10, iii, condition E5']). 
Note that conditions El and E2 hold if R = . 

Let us show that conditions El and E2 hold for "almost all" functions R satisfying 
conditions Rl, R2. More precisely, consider finitely many coupling functions Ri,...,Rn 
satisfying conditions Rl and R2' and take their linear combinations 

N 

Rc(x) = Y / C s R s (x), C=(C 1 ,...,C N )eJR N . 

s=l 

For Rc(x) , conditions Rl and R2' hold if ||C|| < e with a sufficiently small e > 0. Let 
Mi := {C E B £ : condition El holds for R c (x)} and M 2 := {C G B £ : condition E2 holds 
for R c (x)}, where B £ := {C G JR N : ||C|| < e} . In Appendix B, we prove the following 
lemma. 

Lemma 2.11 The sets Mi and M 2 are dense in some ball B £ for a sufficiently small 
e > . 

2.2 Random solution. Convergence to equilibrium 

Let (Q, E, P) be a probability space with expectation E and let B(£ ) denote the Borel 
a -algebra in S . Assume that Yq = Yq(u,p) (see (1.7)) is a measurable random function 
with values in (£, B{£)) . In other words, the map (00, p) 1— > Yq(u,p) is a measurable map 
Q, x IP -> IR 2+2n with respect to the (completed) a -algebra £ x B(JP) and B(M 2+2n ) . 
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Then Y(t) = W(t)Y is also a measurable random function with values in (£,B(£)) owing 
to Proposition 2.7. Denote by fi (dY ) the Borel probability measure in £ giving the 
distribution of Y . Without loss of generality, we can assume (f2, E, P) = (£, 13(E) , /Jo) and 
Y (u,p) = u(p) for p,o(du) x dp-almost all point (uj,p) G £ x IP . 

Definition 2.12 The measure /i t is a Borel probability measure in £ giving the distribution 
of Y(t) , 

m{B) = fM>(W(-t)B), VBeB(£), teM. 

Our main objective is to prove the weak convergence of the measures \i t in the Frechet 
spaces £ s ^ for each s < , (3 < a < —d/2 , 

fJ>t — ? /"oo as t — > oo, (2.8) 
where /ioo is a limit measure on £ = . This is equivalent to the convergence 

| f(Y)nt(dY) - | /(F)/i 00 (rfF) as t - oo 

for any bounded continuous functional f(Y) on £ s ^ . 

Let V = [D F © L> L ] 2 with £> F = C£°(IR d ) , and let D L be the set of vector sequences 
u(k) G IR™ , k G 7L d , such that «(&) = for k G 2 d outside a finite set. For a probability 
measure /i on £ , denote by /t the characteristic functional (Fourier transform) 

jl(Z) = J e i{Y ' z) /2(dY), Z eV. 
Here (•, •) stands for the inner product in L 2 (W) <g> IR^ with different N — 1, 2, . . . , 

(y,z> := E(n^), F = (r°,r 1 ), z = (z°,z 1 ), 

j=0 

(y*, Z { ) := jY\p)Z\p)dp^ J ^(x)C(x)dx+Y.n\k) X l (k), 
P Rd feez d 

where y* = (ip l ,u l ) , = (£\ x*) • A measure /i is said to be Gaussian (with zero expec- 
tation) if its characteristic functional has the form 

jl(Z)=exp{-^Q(Z,Z)}, ZeV, 
where Q is a real nonnegative quadratic form in V . 

Definition 2.13 The correlation functions of the measure /i t , t G IR, are defined by 

Q i t j (p,p') = E(Y i (p,t)^Y j (p',t)), ij = 0,1, P,P eJP, (2.9) 

where E stands for the integral with respect to the measure fio(dY) and the convergence of 
the integral in (2.9) is understood in the sense of distributions, namely, 

(Qy(p,p'),Z 1 (p)®Z 2 (p')) := E{Y i (p,t),Z 1 (p))(Y^(p',t),Z 2 (p')), Z X ,Z 2 G D F @D L . (2.10) 
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2.3 Mixing condition 



Let 0(r) be the set of all pairs of open subsets A, B C IP such that the distance p(A, B) 
is not less than r , and let cr(A) be the a -algebra in S generated by the linear functionals 
Y i — y (Y, Z) for which Z G T> and supp Z C A . Define the Ibragimov-Linnik mixing 
coefficient of a probability measure /i on £ by the formula (cf. [15, Definition 17.2.2]) 



(p(r) = sup sup 



\fj, (A (IB) - h (A)hq(B)\ 



{A,B)eO(r) Aea(A),B ea(B) 
MB) > 

Definition 2.14 A measure fi satisfies the strong uniform Ibragimov-Linnik mixing con- 
dition if ip(r) — > as r — > oo . 

Below we specify the rate of the decay of (p (see Condition S3). 



2.4 Main theorem 

Assume that the initial measure p satisfies the following properties SO S3: 

50. The measure /x has zero expectation value, EY (p) = , p G IP . 

51. The correlation matrices of p, are invariant with respect to translations in Z d , i.e., 
Eqn (1.9) holds for a.a. p,p' G P . 

52. The measure /x has a finite mean "energy" density, i.e., Eqns (1.10), (1.11) hold. 

53. The measure /i satisfies the strong uniform Ibragimov-Linnik mixing condition with 

+oo 

J r d - l V l/2 {r)dr < oo. (2.11) 
o 

Introduce the correlation matrix Qoo(p,p') of the limit measure /i^ . It is translation- 
invariant with respect to translations in Z d , i.e., 

Qoo(p,p') = Qoo(p + k,p' + k), k G Z d . (2.12) 

For Z G V , write 

Qoo(Z,Z) := (Q 00 ip,p'),Zip)^Z{ P ')) = (2n)- d J (^(0), Z n (0, •) ® Z u (9 r ))d9, (2.13) 

where goo(#) is the operator- valued function given by the rule 

a (^VPffl^W+^W) ffW-fflW) ) P(f)) (2U) 

ww-ffw n(eme) + ti\e)) m ' (2 ' 14) 

for G K d \C 1f . Here the symbol ^(0) = Op((?o' ? (^) r ) r ')) stands for the integral operator 
with the integral kernel qQ(0,r,r') (see formula (3.7) with t — 0), r,r' eTZ, and -P/(6*) is 
the spectral projection operator introduced in Corollary 2.9. 

Theorem A Let the conditions SO S3, Rl, R2, El and E2 hold. Then the following 
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assertions are valid. 

(i) The convergence (2.8) holds for any s < and f3 < —d/2 . 

(ii) The limit measure fj,^ is Gaussian on £ . 

(Hi) The characteristic functional of fi^ is Gaussian, 

fc 00 (Z) = exp{-^Q 00 (Z,Z)}, ZeV. 

(iv) The measure ji^ is invariant, i.e., [W(£)]*//<x> = fioo , t G IR . 

Assertions (i)-(iii) of Theorem A follow from Propositions 2.15 and 2.16 below. 

Proposition 2.15 The family of measures {/i t , t G H} is weakly compact in £ s,lS with any 
s < and (3 < a < —d/2 , and the following bounds hold: 

sup£|||W^)y |||L < C{a) < oo. (2.15) 
ten 

Proposition 2.16 The convergence (1.21) holds for every Z G V . 

Proposition 2.15 (Proposition 2.16) provides the existence (the uniqueness) of the limit 
measure //oo • They are proved in Sections 4 and 7, respectively. 

Theorem A (iv) follows from (2.8) because the group W(t) is continuous with respect 
to £ by Proposition 2.7 (ii). 



3 Correlation matrices 

To prove the compactness of the family of measures {fj, t } , we introduce auxiliary notations 
and prove necessary bounds for initial correlation matrices. Since Y l (p, t) = (ip l (x, t),u l (k, t)) , 
we can rewrite formula (2.9) as follows: 

QiW) = E\Y i (p,t)®Yl(j/,t)] = 

Qp\x,x') Qf u \x,k') 
Q^(k,x>) Qf u \k,k>) 



EN)\x, t) ® ip j (x', t)j Ehp^x, t) <g> u j (k', t) 
E\u i {k,t) <g> ip j (x',t)) E[u i {k,t)^u j {k',t)) 



, i,j = 0,1, teJR. (3.1) 



Let us rewrite the correlation matrices Ql J (p,p') by using the condition SI. Note that the 
dynamical group W(t) commutes with the translations in Z d . In this case, condition SI 
implies that 

Q t (k + p,k+p') = Q t (p,p'), t G IR, k G TL d . (3.2) 
Let us introduce the splitting p = k + r , where k G Z d and r G Kf U . In other words, 



x-[x]e Kf, if p = x G TR d , 

o, 

In this notation, (3.2) implies that 



' = i 0, if p = k G 7L d . 
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Using the Zak transform (1.13), introduce the following matrices (cf (2.10)): 

Qi j (9,r,9',r') : = E[Y^(9,r,t) ® Y^(9',r',t)], 9,9' G K d , r, r' G K = T d U 0, (3.4) 

where the convergence of the mathematical expectation is understood in the sense of dis- 
tributions. Namely, write V = [V F @V L ] 2 and V F := C°°(K d x T d ) , V L := [C°°(T d )] n . 
Then 



(Q\ 3 (9, r, 9', /), Zi,(9, r) ® Z^(9', r')) = E(Y^(9, r, t), Z l u (9, r)> (Y^(9\ r>, t), Z> n (9>, r>)) (3.5) 
for Z u = (Zni Z^) eV . Now (3.2) implies that 

Qi j {9, r, 9', r') = {2n) d 5(9 - 9')$ (9, r, r'), 9, 9' G K d , r,r' eK, t G IR, (3.6) 

where 

r, O = e«^>« £ ^(t, r, /) = f ) . (37) 

kez d V Qt \p,y) Qt {V) J 

Recall that IP is the disjoint union IR d U 7L d . For a measurable function Y(p) , write 

J Y(p)dp= J Y(x)dx+ ]T Y{k). (3.8) 
p R d feez d 

Proposition 3.1 Let conditions S0-S3 /io/d. TTien (i) the following bounds hold 

J \Qo(p,p')\ dp < C < oo for any p' G IP, (3.9) 
p 

/ |<2o(p,p')| d P' < C < oo /or any p G IP, (3.10) 
p 

where the constant C does not depend on p,p f G IP . 

(*) V«£,^ j (9,y,yi), V«qf u \9,y), V^(9,y>), (9) are uniformly bounded in 

(9,y,y')eK d xT d xT d , |a|<l-i, \/3\ < 1 - j . 

Proof, (i) By [15, Lemma 17.2.3], conditions SO, S2 and S3 imply 

\Qo(p,p')\ <Cmax{e F ,e L }^ 1/2 (b-p'|), P,p' G P. (3.11) 

Hence, the bounds (3.9) and (3.10) follow from (2.11). 

(ii) Similarly to (3.11), conditions SO, S2 and S3 imply 

K3^( d ^y')\ < E K3^ r (k + y,y')\ <CJ2 <p 1/2 (\k + y-y'\) 

kez d kez d 

< C ]T ^ 1/2 (\k\ -2y/d) < C < oo. 
kez d 

Similar arguments imply the other bounds. ■ 
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Corollary 3.2 By the Shur lemma, Proposition 3.1, (i) implies that the following bound 
holds for any F,GeL 2 := [L 2 (JP,dp)} 2 = [L 2 (IR d ) © [Z 2 (Z d )]«] 2 : 

\(Qo(p,lf),F(p)®Gtf))\<C\\F\\v\\G\\v. 

Corollary 3.3 The quadratic form Qoo(Z, Z) defined in (2.13)-(2.14) is continuous in L 2 . 

Proof. Formulas (2.13) and (2.14) imply 



(Q 00 (p,P , ),Z(p)®Z(p')) = (2n)- d / (q O0 (9,r,r'),Z u (9,r)®Z u (9,r'))d9 



K d \C, 

i oo . 

= (2vr)" d -E J (qo(e, r , r , )+M0,r, r , ),P l (9)Z u (9,r)^P l (9)Z u (9, r '))d9,(3.12) 

l=1 K d \C, 

where f (#,r, r') is the integral kernel of the operator f (9) := 2 n ~,, n ?~nn/n\ 

V ~% (6) n(9)q^(9) J 

Here and below, the symbol (•, •) stands for the inner product in = H Q (Kf) © (D™ , i.e., 
(F,G) = J F~\y)G\y) dy + F 2 " ■ G 2 , F = (F 1 , F 2 ), G = (G 1 , G 2 ) G iff 1 , (3.13) 

K d 

or in H° = [H®] 2 . Consider the terms in the RHS of (3.12). Since 

oo 

E \\Pi z n\\{L 2 (K d x'R,)} 2 ^ C|l^n||[L2 (K d x7e )]2 = C"||Z|| L 2, (3-14) 
i=i 

we obtain 

E J ($(0,ry),Pi(9)zZ(9,r)® Pt(9)Z^(9,r')) d9 < C\\Z\\l, i',f = 0, 1, (3.15) 

l=1 K d \C* 

by Corollary 3.2. Further, consider r l (j . For the terms with fg 1 and r~l° , estimate (3.15) 
holds. We rewrite the term with f® in the form 



E J (r o °°(0, r, r'), fl(0)3£(0, r) ® r')) d9 

l=l K d \C, 

oo „ 

= E 7 (% U (^ r, O, r) ® fi(0)3S(0, d0. (3.16) 

It follows from estimates (3.14) and (9.7) and Corollary 3.2 that the RHS of (3.16) is esti- 
mated from above by C||Z°|| 2 : 2(p) . Finally, consider the term with fg 1 and represent it in 
the form 



oo - _ 

E y (?l\0, r , r >),P l (9)Z 1 u (9, r )®P l (9)P n (9, r >))d9 

l=1 K d \C t 

E y («g°(0, r, r'), n(C)fi(e)Zi(fl, r) ® $1(9)^(9)^(9, r')) d9. (3.17) 



l ~ 1 K d \C 
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By the bounds (3.14) and (9.6), the RHS of (3.17) is bounded by 

c J (||?f>, •, Oilmen-)]" + \\4° u \o, Olltfip?) + IIC V >, Ollflip?) 

M^m)\\Zu(e,-)\\%od0. (3.18) 

In turn, (3.18) is estimated by C||Z 1 ||| 2( -p) by Proposition 3.1, (ii). ■ 

Remark 3.4 The operator qo{6) in H° is nonnegative and self-adjoint. Indeed, for any 
function Z G V , we have 



(q (6), Z u (9, •) ® Z u (6, •)) dO = (2n) d E\(Y, Z) | 2 > 0. 



K d 



Hence, (q (8), Z n (0, •) ® Z n (6, •)) > , G . 

4 Compactness of measures ji t 

Proposition 2.15 follows from the bound (2.15) by the Prokhorov Theorem [25, Lemma II. 3.1] 
by using the method of [25, Theorem XII. 5. 2] because the embedding £°' a C £ s,/3 is compact 
if s < and a > (3 . 

Lemma 4.1 Let conditions S0-S3 hold. Then the bounds (2.15) hold for a < —d/2. 
Proof. Step (i). By condition SI, 

E\\Y(t)\\i, a = s[/(i+N 2 r(i^,oi 2 +iv^,oi 2 +k(x,oi 2 )^ 



|2 



+ j:(i+\k\ 2 r(\u(k,t)\ 2 +\v(k,t)\'< 

kez d 

< C{a,d)e(t), (4.1) 

where 

e(t) :=E[ J (|^(t/,t)| 2 + |V^(i/,t)| 2 + |7r(t/,t)| 2 )^+| M (0,t)| 2 + |t;(0,t)| i 

Kf 

Denote by C t the correlation operator of the random function 

Ti r d\ /tn Tod /r\ rjO/ T^d\ 



y, t) \ yeKf , u(0, t), 7r(y, t) \ yeK{ , v(0, t)) G S 1 := H\Kf) © lR d © H^Kf) © lR d . 

Then e(t) is equal to the trace of the operator C t . Note that C t = Op(q t (0,r,r')^ is an 
integral operator with the integral kernel q t (0,r,r') (see (3.7)), 

9t (0, r, r') = (2vr)- d J e- l(r - r ' )e q t (9, r, r') d0, r, r' G ^ := K{ U {0}. (4.2) 
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Denote by qt(9) := Opyqt(9,r,r')j the integral operator with the integral kernel qt(9,r,r') . 
In this case, (4.2) implies that 



e(t) =tv £l C t = (2n)- d J tr 

K d 



£1 



dO. 



(4.3) 



Step (ii). Introduce the operator T defined by Tipu(0,y) := (V y ipn(9,y),ipn(9,y)) and the 
operator T ex given by 



Then, since qt(9) > and the operator e %rd : (i/)(y),u,ir(y),v) — > (e iy0 ip, u, e iyd 7i, v ) in £ 1 is 
bounded uniformly with respect to 9 G K d , we have 



tr 



e- irB q t {6)e ir ' e < Ctr £l q t {9) = Ctr H o ^^(0)^ 



(4.4) 



Let us now estimate the trace tr H o 
operator Q ex on the space H° , 



^exqt(0)T* ex . Introduce the matrix-valued self-adjoint 



( n(9) o 
V 1 



where I stands for the identity operator on H® . Note that Q ex q t (9)Q ex > (recall that 
Q ex is a self-adjoint operator). Further, B := (T ex Q~ x ) is a bounded operator on H° since 



IT 1 : H°(Kf) © (D n = Hi — > H\ , and ^ ^ J : if* -> 77° . Therefore, 

tr H o[r e ^ t (^r:J = tr H o[77fi e ^ t (0)fi e:r 77*] < Ctr H o[Qe^t(^)^], 



(4.5) 



by [23, Theorem 1.6]. Let us now estimate the trace of the operator Q ex q t (9)Q ex . We first 
use the formula £l ex G(9,t) = U(9,t)Q, ex , where G(9,t) is defined in (11.1) and 



17(0, f) := 



cos fit sin Qt 
- sin fit cos fit 



Hence, by (11.4) we have 

n ex q t {9)Q ex = n ex G(9,t)q (9)G*(9,t)n ex = U(9,t)n ex q (9)Vl ex U*(9,t). 
Since U (t, 9) is a unitary operator on H° , 

tr H o[£Wt(0R*] = ten°[U(9,t)n ex q (9)n ex U*(9,t)] = tr u0 [Q ex q (9)n ex ] 
by [19, Theorem VI.18, (c)]. Finally, it is follows from (4.3)-(4.6) that 

supe(t)<Ci / tv-Ho[^exqo{9)VL ex }d9. 
teR J 



(4.6) 



(4.7) 
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Step (in). Let us now prove that the RHS of (4.7) is finite. We use the representation 

n ex q (9)n* ex = (nexr-^r^^r^^r- 1 )*, 

where T^ 1 stands for the left inverse operator of T ex . On the other hand, Q, ex T~ x is a 
bounded operator in H° since Q ex (9) (T^, 1 , resp.) is (a finite - dimensional perturbation 
of) a pseudodifferential operator of order 1 (—1, resp.) on Kf . Moreover, VL ex {9) is 
uniformly bounded in 9 G K d . Hence, 

tr H o[^o(0R*] < Ctr H o[r e ^ o (0)r:J = CtT £l [q (e)}. (4.8) 

Finally, by inequalities (4.7) and (4.8) and by condition S2, we obtain 

supe(t) < C [ tv £l qo(0)de 

ten J 

K d 

< c&lj (\Mv)\ 2 + IWoO/)l 2 + Mv)\ 2 ) dy + Mo)| 2 + |^ (o)| 2 ; 

Kf 

< d(e F + e L ) < oo. (4.9) 
Now the bound (2.15) follows from (4.1) and (4.9). ■ 

5 "Cutting out" the critical spectrum 

Definition 5.1 i) Introduce the critical set C := C* U ( Ufe Ckj (see El ). 
ii) Introduce the set V° C V given by 

PiZ n (9, •) = for \/l>N,9e K d , 
Zn(6, •) = in a neighborhood of a set CU dK d . 



V° = U N V N , V N := {Z ev 



(5.1) 



Lemma 5.2 Let lim Q t (Z,Z) = Q O0 (Z,Z) for any Z e T>° . Then the convergence holds 
for any Z G V . 

Proof. First, Definition 2.13 implies that 

Q t (Z,Z) := E\(Y(-,t),Z)\ 2 = (Q t (p,p'),Z(p)®Z(p r )), ZeV. (5.2) 
Therefore, by (7.2) we have Q t (Z,Z) = Q (Z(-,t), Z(-,t)) , and hence 

swp\Q t (Z,Z)\ <Csup\\Z(;t)\\h (5-3) 
ten teR 

by Corollary 3.2. By the Parseval identity and by the bound (7.5), we obtain 

\\Z(;t)\\l = C{d) J \\e AT ^Z^{9 r )\\ 2 H>Hl d9<C J \\Z U (9, ■)\\ 2 H o (BHl d9 
K d K d 

= C\\Z\\l (5.4) 

uniformly with respect to t. Here C := L 2 (JR d ) © [l 2 (H d )] n © H\JR d ) © [Z 2 (Z d )] n . Further, 
by Lemma 2.8, for any Z G T> , we can find a Z N G T>^ such that \\Z — Z N \\c — > as 
iV -> oo. Namely, Zg(9) = Ei<n Pi(0)Z n (0) if Z u {9) = for 9 in a neighborhood of 
C U dK d . Finally, the set of such functions Z is dense in C . Then Lemma 5.2 follows from 
(5.3), (5.4), and Corollary 3.3. ■ 
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Lemma 5.3 The convergence (1.21) holds for any ZeD if it holds for Z G T>° . 



Proof. This follows immediately from Lemma 5.2 by the Cauchy-Schwartz inequality: 

\jl t {Z') - h{Z")\ = I / {^ Z>) -eW'%{dY)\ < j \^-z») - l^dY) 



< J | (Y, Z' - Z") \UdY) <sjj | (Y, Z> - Z")\^ t (dY) 
= yjQ t (Z' - Z", Z< - Z") < C\\Z' - Z"\\ c . ■ 

6 Convergence of the covariance 

Proposition 6.1 Let conditions E1-E2, R1-R3 and S0-S3 hold. Then, for any Z G V , 

Q t {Z,Z)^Q O0 {Z,Z), t^oo. (6.1) 

Proof. By Lemma 5.2, it suffices to prove the convergence (6.1) for Z <EV° only. If Z G T>° , 
then Z G V N for some N . Let us apply the Zak transform to the matrix Qt(p,p') , 

(Q t (p,p'),Z(p)®Z(p')) = (27r)- 2d (Q t (e,9',r,r'),Zn(e,r)0Z n (e f ,r')). (6.2) 

Further, by Lemma 2.8, we can choose some smooth branches of the functions Fi(9,r) and 
ui(9) to apply the stationary phase arguments, which requires some smoothness with respect 
to 9. Denote by supp Z n the closure of the set {9 G K d : Z u {9,y) ^ Oj 6 T?} . Since 
supp Zn fl (C U dK d ) = , we can apply Lemma 2.8. Namely, for any point 6 G supp Zn , 
there is a neighborhood (9(0) C K d \ {C U dK d ) with the corresponding properties. Hence, 
supp Zn C u£f =1 (9(0 m ) , where m G supp Zn . Therefore, there is a finite partition of unity 

M 

E 9m(0) = 1, 61 e supp Z n , (6.3) 

m=l 

where g m are nonnegative functions of C^°(K d ) and supp g m C (9(0 m ) . Further, using 
Definition 5.1, ii) and the partition (6.3), represent the RHS of (6.2) as 

M N 

(Q t (p,p'),Z(p) ® Z(p')) = (2n)- d £ E (9m(0)ruit,9),M9) ®A,(0)>, (6.4) 

m=l M'=l 

using formulas (3.6) and (11.6). Here A,(0) = (F,(0, •), Z n (0, •)) , is the 2x2 matrix 

ru>(t,9) := ^E{cos( Wl (0)±^(0))* (^(0) T Q(%hK^)^)) 

+ sinM0)±w,,(0))t (c l {9)p w {9)±p lv {9)CT{9))}, (6.5) 
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where 



(6.6) 



Lur\6) \ / -utf) \ 

J' Cl ^ } -{ur\9) J' 

pgW^^Kfl,-),^)^)^,-)), «eO(e m ), l,V = l,2,..., i,j = 0,1, (6.7) 



and (•,•) stands for the inner product in H® = H°(Tf) © (D™ (see (3.13)) or in H° = 
[H®] 2 . By Lemma 2.8, the eigenvalues ui{9) and the eigenfunctions F t (9,r) are real-analytic 
functions in 9 e supp </ m for every m: we do not mark the functions by the index m to 
simplify the notation. 

Lemma 6.2 Let conditions S0-S3 hold. Then p%{9) e L x (e)(e m )) ; i,j = 0,1, 1,1' = 
1,2,... for each m = 1, . . . , M . 

Proof. Since {Fi{9, •)} is an orthonormal basis, by the Cauchy-Schwartz inequality, we 
have 



m9)\d9 



{f 1 {9,t),^{9,t,t')F v {9,t')) 



d9 



< C J \ P %(9)\ 2 d9< 
o(e m ) o(e m ) o(e m ) 

< J d9 J dr J \qi j {9,r,r')\ 2 dr'. ■ 
o(e m ) 7i 

Further, let us study the terms in (6.4), which are oscillatory integrals with respect to the 
variable 9 . The identities uji(9)+uji>(9) = const + or uji(9)—uji>(9) = const„ with const± ^ 
are impossible by condition E2. Moreover, the oscillatory integrals with u)i(9) ± ui>(9) ^ 
const vanish as t — > oo . Hence, only the integrals with u>i(9) — 0Ji>{9) = contribute 
to the limit because the relation uoi(9) + oj v {9) = would imply the relation uoi(9) = 
ui/(9) = 0, which is impossible by E2. Let us index the eigenvalues u>i(9) as in (2.5). Then 

cos(cj;(#) —ui'(9))t = 1 for 1,1' E (r a -i,r a ], a = 1,2, Hence, for 1,1' e (rv_i, tv] , we 

have 



rw{t,9) 



l(piv(e) + Q(9) Plv (9)Cj(9)) 
+icos2^(^)t -C l {9) Plll {9)C^{9)) 
+ isin2c^)t (Q(0W(^)+PH'(^)Cr(^))- 



(6.8) 



Therefore, 



(Q t (p,p'),Z(p) ® Z(p')} = (2n)- d Y: E fg m (e)(M tt .(e) + ...,M0)®M0))M, (6-9) 
where M w (9) = (M$(0))J =0 , /, /' = 1, 2, . . . , is the matrix with the continuous entries 



M#(0)= XH ,-(F,(0,r), 



5o 



(0, r, r') + Q(9)q (9, r, r')Cf{9)] 13 F v {9, r')); (6.10) 
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here the symbol xw i s given by (see (2.5)) 

J 1 if 1,1' G {r a - U Ta\, a = 1,2, r := 0, , 
Xlv - \ otherwise. ((U1) 

Further, for 9 G suppg m C 0(0) (see Lemma 2.8), we write 

S(^ ) 0=E J? i( 9 - r W(WO. = 0,1. (6.12) 

M'=l 

The local representation (6.12) can be expressed globally in the form (2.14). Hence, 
(Q t (p,p'),Z(p)®Z(p')) = (271)-^ 1 9m(0)(qoc(0,r, r'),Z u (9,r) ®Z n (9,r')) d9 + . . .(6-13) 

where the symbol " . . . " stands for the oscillatory integrals which contain cos(u;i(9)±uJi>(9))t 
and sm(ui(6)±uji'(6))t with oji(9)±uo v (9) ^ const. The oscillatory integrals converge to zero 
by the Lebesgue-Riemann Theorem because the integrands in are summable, and we 
have V(uui(9)±ui'(9)) = on the set of Lebesgue measure zero only. The summability follows 
from Lemma 6.2 because the functions Ai(9) are smooth. The zero-measure condition 
follows as in (2.4) since uoi(9) ± u)i>(9) ^ const. This completes the proof of Proposition 6.1. 



7 Bernstein's argument 

7.1 Oscillatory representation and stationary phase method 

To prove (1.21), we evaluate (Y(-,t),Z) by duality arguments. Namely, introduce the dual 
space £' : = if" 1 - a (IR. <i ) © L' a © H°^ a (JR d ) © L~ a with finite norm 

(III^IIIU) 2 == U\\-i,- a + Iklll-a + IMI-a + \M-a- 
For t G 1R, introduce the "formal adjoint" operator W'(t) , 

(W(t)Y, Z) := (Y, W'(t)Z), Y eS, Z G £', (7.1) 

where (•, •) stands for the inner product in L 2 (TR d ) © [Z 2 (Z d )] n © L 2 (IR d ) © [Z 2 (Z d )] n . 
Write Z(-,t) = W'(t)Z . Then formula (7.1) can be rewritten as 

(Y(t),Z) = (Y ,Z(-,t)), ten. (7.2) 

The adjoint group W'(t) admits a convenient description. 

Lemma 7.1 The action of the group W'(t) coincides with the action of W(t) up to the 
order of components. Namely, W'(t) = exp(^4 T t) , where A is the generator of the group 
W(t) . 
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Proof. Differentiating (7.1) with respect to t for Y,ZeV, we obtain 

(Y,W>(t)Z) = (W(t)Y,Z). (7.3) 

The group W(t) has the generator A (see (1.8)). The generator of W'(t) is the conjugate 
operator 

Hence, A' = A T . U 
Corollary 7.2 The following uniform bound holds: 

||e^)%(Mlk°etfJ ^ C\\Zn(9,-)\\ H » mb Z u (9,-) G H° ® Hi (7.5) 

which can be proved similarly to (9.3). 

Applying Lemma 7.1, we can rewrite Z(t) = W'(t)Z as the Zak transform, i.e., Zji(9, r, t) 
exp (A T {9)t) Zu(0, r). Recall that we can restrict ourselves to elements Z G T>n with a fixed 
index N . Using the partition of unity (6.3), we obtain 

M N 



Z(k + r,t) = (27r)- d ££ [ g m (0)e- i{k+r)e Gf(e,t)F l (e,r)A l (e)de 

m=1 1=1 K d 

N 

= £ £ / e^^W^^^^)*. ZEV N . (7.6) 



Here A t (9) = (F t (9, ■), Z u (9, ■)) , 



Gifrt) := 



/ ,M, sinc^(#)t 

cos ui(9)t — — I 

V ; | , $ G supp g rn , (7.7) 

y — lui(9) sin LUi(9)t cos uoi{9)t 

and c^(#) and af(#) are real-analytic functions in the interior of the set suppg m for every 
m . 

Let us derive formula (1.21) by analyzing the propagation of the solution Z(k + r, t) of 
the form (7.6) in diverse directions k = vt with v G IR d and for r G 1Z . To this end, we 
apply the stationary phase method to the oscillatory integral (7.6) along the rays k = vt , 
t > . Then the phase becomes (Ov ±u>i(9))t , and its stationary points are the solutions of 
the equations v = ^fVuJi(9) . 

Note that Zu(9,r) = at the points (0,r) G K d @1Z with degenerate Hessian Di{9) 
(see El). Therefore, the stationary phase method leads to the following two different types 
of asymptotic behavior of Z(vt,t) as t — > oo . 

I. Let the velocity v be inside the light cone, v = ±Vo^(0) , where 9 G 0(6) \ C . Then 

Z(vt,t) = 0(r d/2 ). (7.8) 

II. Let the velocity v be outside the light cone, v ^ ±Vu>i(9) , where 9 G 0(6) \ C, 
l = l,...,N. Then 

Z(ut,t) = 0(t' k ), Vk>0. (7.9) 
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Lemma 7.3 The following bounds hold for any fixed Z G V° : 

i) sup\Z(p,t)\ < c r d ' 2 . (7.10) 

ii) For any k > , there exist numbers Ct,j > such that 

\Z(p,t)\<C k (l + \p\ + \t\)- k , b|> 7 t (7.11) 

Proof. Consider Z(k + r, t) along each ray k = vt with an arbitrary v G H d . Substituting 
the related expressions into (7.6), we obtain 

N 

Z(vt + r, t) = E / e-^ 6v±UJlW)t e~ t6r a^ (9)^(9, r)A t (9) d9, Z G V N . (7.12) 

m,± 1=1 J d 

This is a sum of oscillatory integrals with phase functions of the form (f)j 1 {9) = 9v ± ui(9) 
and with amplitudes af{9) that are real-analytic functions of 9 in the interiors of the sets 
supp g m . Since uoi{9) is real- analytic, each function <pf has at most finitely many stationary 
points 9 G suppg m (solutions of the equation v = ^fVu>i(9) ). The stationary points are 
nondegenerate for 9 e suppg m by Definition 5.1 and by El since 




(7.13) 



At last, Zu(9,r) is smooth because Z E V . Therefore, we have Z{vt + r,t) = 0(t~ d ^ 2 ) 
according to the standard stationary phase method of [12, 20]. This implies the bounds 

(7.10) in each cone \k\ < ct with any finite c. 

Further, write v := max m max^^jv max |Vu^(0)|. Then, for |t>| > v , there are no 

6»esupp g m 

stationary points in supp Zn . Hence, integration by parts (as in [20]) yields Z{v t + r,t) = 
0(t~ k ) for any k > . On the other hand, the integration by parts in (7.6) implies a similar 
bound, Z(p,t) = 0(^(t/\p\) l ^j for any / > 0. Therefore, relation (7.11) follows with any 
7 > v . This shows that the bounds (7.10) hold everywhere. ■ 

7.2 "Room-corridor" partition 

The remaining constructions in the proof of (1.21) are similar to [5, 8]. However, the proofs 
are not identical, since here we consider a non-translation-invariant case and a coupled 
system. 

Introduce a "room-corridor" partition of the ball {p G IP : \p\ < jt} with 7 taken from 

(7.11) . For t > 0, choose A t and p t G IN . Asymptotic relations between t, A t and p t 
are specified below. Set h t — A t + p t and 

a? = jht, V = a? + A t , j G Z, N t = [{lt)/h t \. (7.14) 

The slabs R\ = {p G IP : \p\ < N t h t , a J < Pd < V} are referred to as "rooms", C{ = {p G 
P : \p\ < N t h u V < p d < a j+1 } as "corridors", and L t = {p G P : \p\ > N t h t } as "tails". 
Here p = (pi, ■ ■ ■ ,Pd) , A t is the width of a room, and p t is that of a corridor. Denote by 
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Xt the indicator of the room R 3 t , by £j the indicator of the corridor C\ , and by r) t the 
indicator of the tail L t . In this case, 

£ t bd(p) + e t (p)] + mi?) = 1, pgp, (7.15) 

N t -1 

where the sum £)t stands for J] . Hence, we obtain the following Bernstein's type 

j=-N t 

representation: 

(Yo, Z(-,t)) = 52 t [{Y ,jdZ(.,t)) + (Y ,e t Z(;t))] + (y ,ifcZ(.,t)>- (7-16) 
Introduce the random variables r\ , c J t and / t by the formulas 

ri = {Y ,)dZ(;t)), 4 = {Yo,&Z(;t)), l t = {Y , Vt Z(;t)). (7.17) 
Then relation (7.16) becomes 

(Y ,Z(.,t)) = J2 t (rl + c{)+k. (7.18) 

Lemma 7.4 Let S0-S3 /ioW and Z e V° . The following bounds hold for t > 1 : 

E\H\ 2 < C(Z) A t /t, Vj, (7.19) 
E\4\ 2 < C(Z) Pt /t, Vj, (7.20) 
^|^| 2 < C7 fc (Z) t' k , \/k > 0. (7.21) 

Proof. Relation (7.21) follows from (7.11) and Proposition 3.1, (i). We discuss (7.19) 
only, and relation (7.20) can be studed in a similar way. Let us express E\r\ | 2 in terms of 
correlation matrices. Definition (7.17) implies 

E\r\\ 2 = (Qo(p,p'),xi(p)Z(p,t)® X i(y)Z(p',t)). (7.22) 

According to (7.10), Eqn (7.22) yields 

E\H\ 2 < Cr d Jxl(p)\\Qo(p,p)\\dpdp' 

= Cr d J xl( P )dp J \\Q (p,p')\\dp'<CA t /t, (7.23) 

where ||(3o(p,p')ll stands for the norm of the matrix (Qo {p,p')^J ■ Therefore, (7.23) follows 
by Corollary 3.2. ■ 

7.3 Proof of Theorem A 

The remaining part of the proof of the convergence (1.21) uses the Ibragimov-Linnik central 
limit theorem [15] and the bounds (7.19)-(7.21). For details, see [5, Sections 8,9] and [9, 
Sections 9,10]. 
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8 Ergodicity and mixing for the limit measures 

The limit measure fj,^ is invariant by Theorem A, (iv). Let be the integral with respect 

to /loo . 

Theorem 8.1 Let the assumptions of Theorem A hold. Then Wit) is mixing with respect 
to the corresponding limit measure , i.e., for any f,g<E L 2 (^, / u oc ) we have 

lim E^fiWWYMY) = E 00 f(Y)E 00 g{Y). (8.1) 
In particular, the group W(t) is ergodic with respect to the measure /i^ , 

1 T 

lim - / f(W(t)Y)dt = E^fiY) [mod ^). (8.2) 



Proof. Step (i). Since ^ is Gaussian, the proof of (8.1) reduces to that of the convergence 

lim E oc {W(t)Y,Z)(Y,Z 1 ) = 0. (8.3) 

t— +OC 

for any Z, Z 1 G V . It suffices to prove relation (8.3) for Z, Z 1 G T>n ■ However, this follows 
from Corollary 3.3 and Theorem A, (iv). 

Step (ii). Let Z, Z 1 G V N . Applying the Zak transform and the Parseval identity, we 
obtain 

I(t) = E 00 (W(t)Y,Z)(Y,Z 1 ) = E 00 (Y,W'(t)Z)(Y,Z 1 ) 

= {2n)- d {q oo {0, r, r'), G*(9, t)Z n (9, r) ® 1^(9, r')). (8.4) 

Using a finite partition of unity (6.3), and relations (8.4) and (6.12), we see that 

M N . 

7(0 = (2vr)- d E gm{9)G* l {9,t)A l {9)M w {9)A v {9)d9. (8.5) 

m=l 1,1' =1 J 

Here A t (9) = (F,(0 , •) , Z n (9 , •)) and A v {9) = {^(9, ■), Z^(9, ■)) , and Gi(t,9) is defined in 
(7.7). Similarly to (7.6), we have 

Ht) = E E / 9m{9)e ± ^ t af{9)A l {9)M l v{9)A v {9) d9. (8.6) 

m l,l'=l J 

Here all the phase functions ui(9) and the amplitudes af{9) are smooth functions on 
suppg m . Further, the relation Vuji{9) = holds on a set of Lebesgue measure zero only. 
This follows similarly to (2.4) since Vuji{9) ^ const by condition E2. Hence, I(t) — > as 
t — > oo by the Lebesgue- Riemann theorem since the functions M U i{9) are continuous. ■ 

Remark. A similar result for wave equations and for harmonic crystals was proved in [4, 8]. 
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9 Appendix A: Dynamics in the Bloch-Fourier repre- 
sentation 



In this appendix we prove the bound (2.2). We first construct the exponential exp(A(9)t) 
for any chosen 9 G K d = [0, 2ii] d and study its properties. Let us choose 9 G K d and 
X G H 1 := H{ © if° , where H[ = H s (Tf) © (D ra . Introduce the functions exp(A(9)t)X 
for X G H 1 as the solutions X(9,t) to the problem 

/ X(9,t)=A(9)X(9,t), teJR, , , 

\ X(9,0)=X . iy ' iJ 

Proposition 9.1 For any chosen 9 G K d , the Cauchy problem (9.1) admits a unique solu- 
tion X(9,t) G C(JR;U r ) . Moreover, 

X(9,t) = e A{e)t X , (9.2) 

and 

||*(M)IIhi <C||X || H i, (9.3) 
where the constant C does not depend on 9 G K d and t G 1R . 

We prove this proposition in Subsection 9.2. 



9.1 Schrodinger operator 

Let us first construct solutions X{9,i) to problem (9.1) with a chosen parameter 9 G 
K d . Write X(9,t) = (X°(6, t), X\0, t)) , where X°(9,t) = (ip(9,t),u(9,t)) and X 1 (9,t) = 
(<f)(9,t),v(9,t)). By (9.1) and (1.15) we have X\9, t) = X°(9, t) , and X°(9,t) is a solution 
to the following Cauchy problem with a chosen parameter 9 G K d : 



X°(9,t) = -H(9)X°(9,t), teR, 
(x°(9,t),X%9,t))\ t=Q = (X °,X 1 )=X , 



where £1(9) = \JH{9) > . 

Lemma 9.2 For X° G H® , the following bounds hold: 



where the constant C does not depend on 9 G K d . 



(9.4) 



where H(9) is the "Schrodinger operator" (1.16). Hence, formally, 

x°(e,t) = co S n(9)tx° + smn(9)tn- 1 (9)x^, (9.5) 



M9)X \\ H -, < C\\X°\\ H o, (9.6) 
\\n- 1 (9)X°\\ H o < C\\X°\\ H o, (9.7) 
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Proof (i) Formula (1.16) for H(9) implies that 

\\H(6)X°\\ H -,<C\\X°\\ Hl , X°eH±. (9.8) 

where the constant C does not depend on 9 G K d . Hence, 

\\mxY H » = (X°,H(9)X°) < ll^llfljIlW^ll^i < C\\XX Hl - (9.9) 

Since, Vt(9) = SI* (6) , the bound (9.9) implies (9.6). 
(ii) Condition R2 implies that 

IIX^JH^X ^ > (X°,H(6)X°) > k 2 \\XYhI- 

Hence, \\H(9)X°\\ H -i > k 2 ||X°|| h i. Therefore, \\H- 1 (9)X°\\ H i < k' 2 \\X°\\ h -i. In particu- 
lar, (9.7) follows. 1 ■ 

Remark 9.3 Condition R2' implies condition R2. 

Proof. Indeed, for X° = (ip(y),u) G Hi and 9 G K d we have 

(X°,H(9)X°) = 



ip{y)[(iV y + &Y + m 2 My) + Rn{9, y) ■ u(9)ip(y) 



+Rn{e,y)-u{d)Tp{y) dy + cu 2 (9)\u 



^ y + 9My)f + ^ 



cp(y) + —R n (9,y)u(9) 
m 



dy 



+u 2 (9)W^ 2 



> a 



R u (9,y)u(9)\ 2 dy + ^ J My)\ 2 dy 
J \V y ^y)\ 2 dy+(^-Pd2 7 r) J \<p(y)\ 2 dy 



T d 



+H9)\ 2 (u 2 -^J\Rn(9,y)\ 2 dy) 
T d 



for some (5 > and a G (0, /3/(/3 + 1)) . Take j3 < m^j (And) . It remains to prove that 

2 



vl-—\ \Rn{9,y)\ 2 dy>Q. 



T d 



With regard to condition R2', The Parseval equality implies 

J \R n (9,y)\ 2 dy = J | £ e ikB R(k + y)f dy < J | £ \R(k + y)( dy < u 2 m 2 /2. 

rpd rpd fegZ d rpd k^TJ 1 
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9.2 Existence of the Schrodinger group 

Recall that u>i(9) > ( Fi(9, •) ), I — 1,2, ... , are the eigenvalues (orthonormal eigenvectors) 
of the operator Q(9) in H® . Let us prove the existence of solutions to the Cauchy problem 
(9.1). We represent X°(9,t) in the form 

oo 

X°(9,t) = J2Mt)Fi(9,r), teJR, (9.10) 
i=i 

where Ai(t) = Ai(9,t) is the unique solution of the Cauchy problem 

Mt) = -cof(9)Mt), (Mt),Mt))\ t=0 = (A^A^), 
and Ai a = A* a (e) = (F l (6,-),X*(-)), % = 0, 1 . Hence, 

A l (t) = cos^0)tA 0l + S ^^Al (9.11) 

By the energy conservation, this yields 

\Mt)\\ , (0) \Mt)\ 2 _\A l 0l \ 2 u 2 (Q Mli? 



2 

Summing up, for t G IR , we obtain 
\\\X"{9,t)f HO + \{X\9,t), H(0)X°(e,t)) = \\\Xl\\l + \{XlH{9)X°) < C\\X Q \\^ (9.12) 

by (9.8). Hence, the solution (9.10) exists and is unique. 

Further, relation (9.11) implies (9.5). Finally, the solution to problem (9.1) exists; it is 
unique and can be represented by (9.2). The bound (9.3) follows from (9.12) and (2.1). ■ 

Now the exponential exp (^A(9)tj is defined for any chosen value 9 G K d , and this 
exponential is a continuous operator in H 1 . 

9.3 Smoothness of the Schrodinger group 

To complete the proof of Proposition 9.1, we must prove the smoothness of the exponential 
with respect to 9 . This is needed to define the product (1.19) of the exponential and the 
distribution Y on (-) . 

Consider the operators exp (A'(9)tj , t E IR, on H" 1 := (H 1 )* = iff 1 © , where 
A!{9) is the formal adjoint operator to A{9) : 

(X,A'(9)Z) H o = (A(9)X, Z) H o, X,Ze C™(T?) x <D». 

Note that 

A'(9)=A T (9)=(° 1 ~ 7i{9) ). (9.13) 
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Lemma 9.4 For any a > , the following bound holds: 

sup sup ]T IIP^'^XoIIh-i <C(T)||X ||h-i. (9.14) 

|i|<Tee^ d | 7 |< a 

Proof. For ct = , the bound 

||e* w %|| H -i ^ CII^oIIh-!, (9.15) 

follows from the bound (9.3) by duality arguments. Consider the case a — 1 . Introduce the 
function X 7 (t) := V]X(6,t) , where X(6,t) = e^'W*X . Then 

X 7 (t) = ^L'(0)X 7 (t) + [D^'(^)]X(^, i), X 7 (0) = 0. 

Hence, 

t 

X 7 (t) = | e^'W( t - s )[I?^ / (^)]X(^,s)rfs. 
o 

Therefore, by the bound (9.15), 

t 

||X 7 (t)|| H -i < f\\e*WM[V2AWX(6,8)\\ s -id8 



t 

< Cj \\[V]A\e)]X(e,s)\\ H ^ds. (9.16) 



It follows from (9.13) that 

[VeA m ~ [ J ' [P ^ WJ - [X>?S*(6>)] 2^{9)V^{9) ) ' 

Here [P?S(0)] u := [P^ n (0,-)]«, «6C",and [P^*(0)] y^) := J [V] R n (-9 , y)]<p(y) dy . 
Hence, if 9 e K d , then 

^'(^(MIIh- = ^(^(Mll*- < ch^ 1 ^^)!!^ 

< C7||e*W'X ||e-i < C||X || H -i, (9.17) 

by the bound (9.15). Inequalities (9.16) and (9.17) imply the bound (9.14) with a — 1 . For 
a > 1 , the estimate follows by induction. ■ 

9.4 Dual group 

Here we complete the proof of the bound (2.2) by duality arguments. Introduce the dual 
space £' := H~ 1 - a (Wl d ) © L~ a © H°'~ a (lR d ) © of functions Z with finite norm 



i^iio,-j 2: = iHr-i,- Q +iK- Q +Ni^+ iL ' 112 



V 



For Z G 5' , we have Z n (0, •) G H a (K d ; H" 1 ) . 
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Lemma 9.5 Let a be even and let a < —2 . Then 



siivwm%,- a <c(T)\mo,- a - (9- 18 ) 

\t\<T 



Proof. Note first that 



(\\\Z\\\' ,- a ) 2 - E / mZuWr^M (9.19) 

h\<-a K d 

for Ze8'. Indeed, 

ll«l|ia = E (k)- 2 >(k)\ 2 = C[\(1- Ae)- a/2 u(9)\ 2 d9, 
kez d T d 

li,-a = \\(x)- a ^{x)\\h - HA-V)-^)!!!, = c||(i + iei 2 )- 1/2 (i-A,)-^(Olli 2 

E / (1 + \2nm + ei 2 )- 1 / 2 ^ - A e )- a/2 ij(2nm + 9)\ 2 d9 



meZ d Kd 



E / (l + |m| 2 )- 1 / 2 |(l-A e )- Q / 2 ^(27rm + ^)| 2 ^ 



meZ d Kd 

~ /||(l-A e )-^n(^-)ll^)^- (9-20) 

Hence, by Lemma 9.4 and by (9.19), 

(\\\W'(t)Z\\\'^ a ) 2 ~ E /||2?j(^Z n (fl, 0)ll^-i d» 



l7l<-a^d 



<C(f) E /ll^n(fl,0||^-ide~C7(0(lZK > _ a ) 2 . ■ 
Corollary 9.6 T/ie bound (2.2) follows from (9.18) by the duality considerations. 

10 Appendix B: Crossing points 
10.1 Proof of Lemmas 2.8 and 2.10 

Let us prove Lemma 2.8. For any chosen 9 G IR d , the Schrodinger operator H(9) admits 
the spectral resolution 

oo 

i=i 

where < Xi(0) < M{9) < and P/(0) are one-dimensional orthogonal projectors in 
. Further, let us take an arbitrary point 9 G IR d and a number A G (Am (6), Am+i(©)) , 
where M > N and A M (6) < A M+1 (6) . Then A ^ A,(0) for G (9(6) if (9(6) is a 
sufficiently small neighborhood of 6 . Write 

n A (o) = E \(0)Pi(9) 

A ; (e)<A 

P A (0) = E Pi(0) 

A,(e)<A 



G (9(0). 
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Further, let us choose a contour T\ (in the complex plane (D ) surrounding the interval 
(0, A) such that AgTa. In this case, by the Cauchy theorem, 



H A (0) = 



P A (0) 



\d\ 



r A 



H(0) - A 



dX 



r A 



7i(0) - X 



G (9(6). 



Finally, by (1.16)-(1.18) and by condition Rl, the mapping 71.(0) is an analytic operator- 
valued function in G O c (Q) , where O c (Q) is a complex neighborhood of (9(6) . Therefore, 
the same integrals converge for G O c (Q) , and the functions P A (0) and 7i A (0) are an- 
alytic in a smaller neighborhood 0' C (Q) . Reducing 0' C (Q) again, we can choose a basis 
e\(0), . . . , eu(0) in the space R A (0) := P A (0)H® , where the functions ei(0) depend analyt- 
ically on G O c (Q) . For example, it suffices to choose an arbitrary basis ei(6), . . . , e M (Q) 
and set e t (0) = P A (0)ei(Q) . The operator 7t(0) on the invariant space R A (0) can be 
identified with the corresponding matrix 



n A (o) = (n kl ( 



k,l=0,...,M : 



;io.i) 



which depends analytically on G O c (&) . Therefore, the eigenvalues X±(0), . . . , Xm(0) 
and eigenvectors Fi(0), . . . ,Fm(0) of this matrix can be chosen as real- analytic functions 
of G O r (Q) \ C A , where O r (Q) := O c (Q) n H d and C A is a subset of IR d of Lebesgue 
measure zero. This can be proved by using the methods of [8, Appendix]. It remains to 
pass to the limit as A — > oo and define C* := U^°C A . Finally, u>i(0) := \J\i(0) . After this, 
relations (2.5) and (2.6) follow as in [8, Appendix]. ■ 
Lemma 2.10 can be proved in a similar way. 



10.2 Proof of Lemma 2.11 

First let us show that conditions El, E2 hold for Rc(x) = corresponding to C± — . . . — 
C N = . Indeed, in this case, relation (1.16) becomes 




H(0) := 

Therefore, ui(0) are equal to either u*(Q) or yj (2irk + 0) 2 + wiq , k G 7L d . Namely, u*(6) 
corresponds to the eigenvectors (0, u) with an arbitrary u G IR™ . The square root corre- 
sponds to the eigenvectors F k (0,y) = (e~ 2mk ' y ,0) with k G Z d . It can be readily be seen 
that conditions El and E2 hold in this case. 

Further, choose an arbitrary / = 1,2,..., a point G IR d \ C* and a bound A G 
(Xm(0),Xm+i(0)) as above (with M > I). The function Rc(x) and the corresponding 
operator 71q(0) depend analytically on (0, C) G (C d x (C N . Moreover, Rc(x) satisfies 
conditions Rl and R2' for C G B e with a sufficiently small e > 0. Therefore, as in the 
proof of Lemma 2.8, the corresponding eigenvalues ui(0,C) , I — 1, . . . , M , are also analytic 
functions of (0, C) in the domain A4i(Q) — O c \C , where O c is a complex neighborhood of 
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O r (Q) x B £ , and C is a proper analytic subset of O c . Hence, the corresponding determinant 
Di(9,C) is an analytic function of Mi(Q) . Further, Mi(Q) is an open connected set since 
C is a proper analytic subset. Therefore, Di(9,C) ^ on Mi(Q) since Di(9, 0) ^ 0, 
9 G !R d . Further, introduce the set 

M ll = {CeB e :D l (6,C)£0}. 

The set B £ \ Mu cannot contain any open ball, since otherwise Di(9,C) = 0. Hence, M\\ 
is an open dense set in B £ . It remains to note that Mi = fl/Mi/ is thus a dense subset of 
B £ . For M 2 , the proof is similar. ■ 



11 Appendix C: Covariance in the spectral represen- 
tation 

Introduce the matrix-valued operator 

y — \ 1(9) smi 1(9 )t cosil(9)t J 
Note that we can represent the matrix G{9, t) in the form 

G(9, t) = cos Q(9)t I + sin Q(9)t C(9), (11.2) 
where / stands for the unit matrix, and 

Q-\9) \ 



C(9) 



-n(e) o 



In this case, the solution of (1.14) has the form Yn(9,r, t) = G(9,t)Y n(9,r) , r e 1Z . Using 
(11.2) and (3.4), we obtain 



Q t (9,r,9',r') = E[Y u (9,r,t) ®Y n (9',r',t)} 

= cos Q(9)t Q (9, r, 9', r') cos fl(9')t 

+ sintt{9)t C{9)Q (9,r,9',r')C T (9') smQ{9')t 

+ cosQ(9)t Q {9,r,9',r')C T {9') smQ(9')t 

+ sin Q(9)t C(9)Q (9, r, 9', r') cos Q(9')t. (11.3) 

By (3.6), we see that 

q t (9) = G(9,t)q (9)G*(9,t) = cosQ{9)t q (9) cosQ(9)t 

+ cosn(9)t q (9)C T (9) sinn(9)t + sin n(9)t C(9)q (9) co S n(9)t 

+ smQ(9)t C(9)q (9)C T (9) sintt(9)t, (11.4) 

where q\{9) is the integral operator with the kernel qt(9,r,r') defined by (3.7). 

For the simplicity of our manipulations, we assume now that the set of "crossing" points 
is empty, i.e., u>i(9) ^ u>i>(9) for any 1,1' e IN, and the functions u>i(0) and Fi(9,r) are 
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real- analytic. (Otherwise we need a partition of unity (6.3)). Consider the first term in the 
RHS of (11.4) and represent it in the form 

cosn(6)t q (9) cosQ(9)t = F t (9, r) (cosc^(0)i p w {9) ca&u v (6)i)Fv(0, r') 

i,i' 

= ^FK^r)i[cos(^K^)-^K^))t + cos(^K^)+^(^))t]^WWO, (H-5) 
i,i' 1 

where p iv {9) = (v%(0))) j=q = {F l {9, •), ($(e)F v )(6, •))* . =Q (p%(0) are introduced in (6.7)). 
Similarly we can rewrite the remaining three terms in the RHS of (11.4). Finally, 

oo 

$ j (9,r,r') := £ r)rg(f, 0) ® *R0, r'), (11.6) 
M'=i 

where r U t(t,9) = (r'lj,(t,9))} j=0 are defined in (6.5). 
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